Theory of two-magnon Raman scattering in iron pnictides and chalcogenides 
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Although the parent iron-based pnictides and chalcogenides are itinerant antiferromagnets, the use 
of local moment picture to understand their magnetic properties is still widespread. We study mag- 
netic Raman scattering from a local moment perspective for various quantum spin models proposed 
for this new class of superconductors. These models vary greatly in the level of magnetic frustration 
and show a vastly different two-magnon Raman response. Light scattering by two-magnon excita- 
tions thus provides a robust and independent measure of the underlying spin interactions. In accord 
with other recent experiments, our results indicate that the amount of magnetic frustration in these 
systems may be small. 
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Since their discovery in 2008 [1, 2], iron-based su- 
perconductors have led to extensive studies in the con- 
densed matter field. Like the high-Tc cuprate supercon- 
ductors, the iron pnictides and chalcogenides are quasi- 
two-dimensional materials with a layered structure. Both 
the copper and iron based superconductors exhibit anti- 
ferromagnetic (AF) order in the parent phases, and su- 
perconductivity emerges by suppressing magnetism upon 
doping. However, there are also clear differences. At sto- 
ichiometry, the cuprates are Mott insulators and can be 
described by an effective single-band model. On the other 
hand, the iron-based pnictide and chalcogenide parent 
compounds are multi-orbital systems and remain metal- 
lic even in the AF state. The description of magnetism 
and its interplay with superconductivity in this new class 
of superconductors remain an active area of research. 

There is growing evidence that weak coupling the- 
ories with varying degrees of sophistication and input 
from density functional theory (DFT) calculations can 
quantitatively describe several experimental findings [3- 
5]. However, many aspects of the magnetic properties in 
these materials are also remarkably captured by a local 
moment perspective. These include the phase diagram 
with orthorhombic distortion and antiferromagnetic or- 
der, temperature dependence of uniform susceptibility, 
and neutron studies of spin-wave dispersion throughout 
the Brillouin zone (BZ). It also has been shown that a 
local moment model can reproduce the essence of mag- 
netism from DFT calculations [6]. There are various 
reasons why a strong coupling perspective is still rele- 
vant. It provides the potential for a unified framework for 
understanding high-Tc superconductivity derived from 
electron-electron interactions. More importantly, Local 
moment models are being widely used in experiments to 
explain magnetic properties in these systems [7-13]. 

In the pnictides, proposed theories for their coUinear 



AF phase based on localized moments include the spa- 
tially anisotropic Jia-Jib-J2 model [14, 15], where the 
coupling is ferromagnetic (FM) in one direction and AF 
in the other, and the strongly frustrated J1-J2 model [16- 
19], where coUinear antiferromagnetism arises via order 
by disorder (see Fig. 1). Similarly, the diagonal double 
stripe AF order in iron chalcogenides can be obtained by 
either invoking strong frustration in a J1-J2-J3 model [20] 
or utilizing a model with strong spatial anisotropy stem- 
ming from orbital order [21]. These Hamiltonians can 
lead to ground states of the same broken symmetry, but 
they vary greatly in the degree of magnetic frustration. 
It is necessary to distinguish these scenarios and narrow 
the possible models for the iron-based superconductors. 

The preceding difficulty in some cases could be resolved 
by inelastic neutron scattering (INS) experiments. In- 
deed, the INS spin wave spectrum of CaFe2As2 [7] has a 
local maximum at momentum (tt, tt) (in the single iron 
BZ), strongly favoring the Jia-Jib-J2 model. However, as 
high energy spin wave signals are strongly damped, in- 
terpretations of the neutron results remain disputed [22] . 
Moreover, there are also cases where different Hamilto- 
nians can lead to virtually indistinguishable spin wave 
spectra. Therefore, further studies based on different 
measurements appear crucial. 

One way to achieve the above goal is by studying mag- 
netic Raman scattering [23]. Light scattering by two- 
magnon flips is dominated by short-range excitations and 
sensitive to details of the exchange couplings. This probe 
was instrumental in the first accurate determination of 
the exchange constants in cuprate superconductors [24]. 

In the present study, following the Fleury-Loudon (FL) 
formalism [25, 26], we investigate the two-magnon Ra- 
man spectra for strongly frustrated models and spatially 
anisotropic models that are unfrustrated. Using exact 
diagonalization (ED) and spin wave theory, we find that 
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FIG. 1. (Color online) Schematics of quantum spin models 
that lead to (a) the (tt, 0) coUinear AF order in iron pnictides, 
and (b) the (■7r/2, 7r/2) diagonal double stripe AF order in iron 
chalcogenides. 



in strongly frustrated models the two-magnon peak ap- 
pears at relatively low energies. Strong frustration im- 
plies abundant low-energy configurations, which thereby 
provides symmetry-dependent pathways for low-energy 
two-magnon flips. In contrast, in unfrustrated models 
the two-magnon peak occurs at a higher energy, closer to 
twice the single-magnon bandwidth. Magnetic Raman 
scattering thus can be a robust and independent mea- 
sure of the underlying spin interactions. 

We study quantum spin models on a square lattice: 



(1) 



where Jij are short range exchange couplings up to the 
third nearest-neighbor (NN). Depending on the interac- 
tion parameters, the model can support ground states of 
different broken symmetries, such as the (tt, 0) coUinear 
AF order or the (7r/2,7r/2) diagonal double stripe order. 
The FL Raman scattering operator O is given by [23] 



i,3 



(2) 



where the coupling strengths rjij are proportional to the 
exchange interactions Jy , dy is the unit vector that 
connects lattice sites i and j, and Gin and Gout are 
the incident and scattered photon polarizations, respec- 
tively. Specifically, we will consider the following light- 
polarization geometries: 

■^{x + y),eout = -^(x-y), 
-i=(£ + y),eo„t = -^(x-hy), (3) 
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where we note that A'^^ transforms as Aig © B2g [23]. 
With a given scattering operator O, the two-magnon Ra- 
man cross-section is given by = — (l/7r)Im[/(aj)], 
where 



1 



Eo + ie-H 



(4) 



Here |^'o) is the ground state with energy Eq. 

Below we perform spin one-half {S = 1/2) calculations 
in systems with a (tt, 0) coUinear AF order, and spin one 
{S — 1) calculations in systems with a (7r/2,7r/2) di- 
agonal double stripe order. These calculations are rele- 
vant since in the pnictides the Fe moment is unexpectedly 
small (ranging from 0.3 to 0.9 hb), while in the chalco- 
genides the moments (~ 2/xb) are much larger [27]. For 
the S' = 1/2 case, we perform ED calculations on square 
lattices of 16, 32, and 36 sites. For the S = 1 case, we 
focus on a 16-site cluster as the Hilbert space is substan- 
tially larger. We further compare the calculations with 
spin wave theory and comment on the effects of 5*. 

We first discuss the case of (tt, 0) AF order. We con- 
sider (i) a frustrated Ji — J2 model with Ji = J2 and (ii) 
the spatially anisotropic Jia-Jib-J2 model, as depicted 
in Fig. 1(a). For the latter we take Jib = — O.lJia 
and J2 — 0.4Jia as reported by INS measurements on 
CaFe2As2 [7]. In this case, even though [ Jia -f Ju) is 
comparable to 2J2, the system remains unfrustrated. 

Figure 2 shows the two-magnon Raman spectra. In the 
36-site ED calculations, the J1-J2 model has a Big peak 
at ~ SJi, while in the A'-^g channel it is at ^ 65' Ji. On 
the other hand, for the Jia-Jib-J2 model the two-magnon 
excitation appears at ~ ISJia in both polarizations. We 
note that while finite-size effects give small corrections 
to the Raman resonance energy, the clear qualitative dif- 
ference in the Big channel between the two models is 
robust, depending weakly on cluster sizes. 

The sharp distinction in the Big two-magnon response 
is related to the difference in the amount of magnetic frus- 
tration. In both models, the single-magnon bandwidth 
is roughly the same [29]; however, their magnon spectral 
weights are rather distinct. In the Jia-Jib-J2 model, the 
magnon density of states (DOS) diverges at the top of 
the single-magnon band, similar to that in the Heisen- 
berg model. On the other hand, in the frustrated J1-J2 
model a large spectral weight is distributed at lower en- 
ergies [28]. In general, a substantial low-energy magnon 
DOS as implied by frustration provides pathways for low- 
energy two-magnon flips [30, 31]. 

The low-energy Big mode of the J1-J2 model can be re- 
lated to the proximity of the system to a disordered phase 
characterized by either a quantum spin liquid or a valence 
bond solid (VBS) . At a transition to a phase lacking mag- 
netic order, one expects singlet modes with appropriate 
quantum numbers to soften to zero energy [32]. Close 
to such a critical point low energy excitations would ap- 
pear. For a quantum spin liquid, low energy singlets may 
appear in all light polarizations, while for a VBS phase 
low energy singlets may appear only in specific channels. 
Our finding of a low-energy excitation in Big and not in 
A'lg implies a nearby VBS phase with columnar dimcrs, 
since this state has the right symmetry [33, 34]. 

To compare the calculations to experiments, we con- 
vert our results into physically measurable units. By 
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FIG. 2. (Color online) The Fleury-Loudon two-magnon Ra- 
man cross-sections calculated with S = 1/2 exact diagonal- 
ization. In the frustrated J1-J2 model, the Big two-magnon 
excitation appears at ~ SJ\ , an energy much lower than twice 
the single-magnon bandwidth. 



fitting to the Jia-Jib-J2 model, an SJia ^ 50 meV 
has been reported by neutron scattering experiments 
on CaFe2As2 [7]. Therefore, the 36-site ED calcula- 
tion indicates a Big two-magnon peak at lAJia ^ 3000 
cm~^. Anisotropic exchange couplings are also found 
in BaFe2As2 [8] and SrFe2As2 [9], with a slightly differ- 
ent energy scale. Even with a material dependence and 
possible corrections due to finite-size calculations, the 
two-magnon peak in an unfrustrated Jia-Jib-J2 model is 
clearly at an energy several times higher than SJia- On 
the other hand, for a J1-J2 model with the same exchange 
energy scale, a Big two-magnon peak occurs around a few 
hundred wavenumbers. As in the cuprates [35], the ab- 
sence of a low-energy Big two-magnon peak in recent ex- 
periments implies that magnetic frustration in the pnic- 
tides might also be small [10, 11]. 

To further understand the Raman cross-sections, we 
discuss results based on spin wave theory. In the AF 
Heisenberg model, linear spin-wave theory indicates a 
Big two-magnon peak located at 8SJ, twice the energy of 
the single-magnon bandwidth. This value is appropriate 
only in the classical S = 00 limit, where magnon-magnon 
interactions can be neglected. In systems with a finite S, 
two magnons are bound to each other locally in space, 
thereby shifting the two-magnon energy to a lower value. 
When interaction effects are treated at the mean-field 
level, the two-magnon Raman response takes the follow- 



Im 



jM/(i-k|/H) 



where 



ing RPA form: 

I{lu) [defined in Eq. (4)] is calculated within the linear 
spin wave approximation, and J is the mean-field inter- 
action characteristic to the problem. In the 5 = 1/2 
AF Heisenberg model, the Big two-magnon energy is 
renormalized to ~ 6SJ. Similarly, in the 5 = 1/2 ED 
calculations for the Jia-Jib-J2 model, the two-magnon 



excitation in the Big channel occurs at ^ 7 SJia, again 
understood by an RPA renormalization of its classical 
Big two-magnon peak at SSJia- 

In frustrated models, however, the correspondence be- 
tween the classical Big two-magnon energy and twice 
the single-magnon bandwidth fails to describe the Ra- 
man spectra. Based on linear spin wave theory, while 
the top of the single magnon band is roughly 65Ji in 
the Ji- J2 model, its classical Big two-magnon energy oc- 
curs at zero energy. This zero energy two-magnon res- 
onance is closely related to the zero mode in the spin 
wave spectrum at momentum (7r,7r), which is known to 
be shifted to a finite energy by quantum fluctuations. 
The spin- wave theory calculation agrees well with the ED 
results shown in Fig. 2, where the Big two-magnon en- 
ergy ~ SJi is much closer to zero, rather than twice the 
single-magnon bandwidth. Moreover, our 5 = 1 16-site 
ED calculation indicates that the Big two-magnon peak 
moves further down to ^ 0.75Ji. The distinction in the 
two-magnon response thus serves as a clear benchmark 
to distinguish the Ji- J2 and the Jia- Jit- J2 models. 

We next focus on situations where the model Hamil- 
tonian [Eq. (1)] has a (7r/2,7r/2) AF ground state as 
observed in iron chalcogenides [36, 37]. Like the case of 
the pnictides, below we discuss two models that have a 
similar single-magnon bandwidth but vary in the amount 
of magnetic frustration. 

In the first case we consider a system which is spatially 
anisotropic but unfrustrated. We use parameters from 
DFT calculations [15] (in units of meV): SJia = —7.6, 
5Jih = -26.5, 5J2a = 46.5, and 5J2b = -34.9 [see Fig. 
1(b)]. In the second frustrated case, we use parameters 
from Fei.osTe INS measurements [12] (in units of meV): 
SJia = -17.5, SJib = -51.0, 5J2a = 21.7, 5J2b = 21.7, 
and 5J3 = 6.8. In the latter case a strong frustration 
is present, as the NN interactions in both crystal- axes 
are FM, and the second NN exchanges are all AF, inde- 
pendent of spin sub-lattice. In both cases, SJ below is 
referred to as the energy scale of the dominant coupling. 

Figure 3 shows the two-magnon Raman spectra from 
16-site 5=1 ED calculations. In the unfrustrated sys- 
tem, there is no low-energy two-magnon excitation. The 
two-magnon Raman spectra start to show spectral fea- 
tures between 6 ~ 85J, and the cross-sections are dom- 
inated by peaks located at ^ 115J and 10.55J respec- 
tively in the Big and A[g polarizations [Fig. 3(a) and 
(b)]. In this unfrustrated case, the two-magnon peak is 
expected to occur at an energy close to twice the single- 
magnon bandwidth [see Fig. 4]. On the other hand, when 
the system is frustrated there is a substantial low-energy 
two-magnon weight between 1 ~ 35J [Fig. 3(c) and (d)] 
These results show a clear difference in the two-magnon 
Raman response between the two models. 

As mentioned previously, the distinction in the two- 
magnon Raman spectra is not caused by a difference in 
the single-magnon bandwidth. As shown in Fig. 4, the 
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FIG. 3. (Color online) Two-magnon Raman spectra for 
models with a (7r/2,7r/2) diagonal double stripe AF order. 
In both cases SJ is the energy of the dominant exchange 
coupling. We note the plots have a different energy scale. 



magnon dispersion from spin wave theory [12, 20] and 
the dynamical structure factor S^^{q,uj) from ED calcu- 
lations both indicate a similar magnon bandwidth of the 
two models. However, in a frustrated model the magnon 
DOS is no longer dominated by zone boundary magnons, 
and an estimate based on twice the magnon bandwidth 
fails to describe the two-magnon Raman profile. 

We last note that recent Raman scattering experiments 
do not suggest the existence of low-energy two-magnon 
excitations in the chalcogenides [13]. This implies that 
weakly frustrated but spatially anisotropic models based 
on orbital order [21] might be relevant. When there is 
an ambiguity or freedom in fitting the neutron scattering 
spin wave spectra to localized model models, additional 
experiments such as Raman scattering can be helpful for 
pinning down the interaction parameters. 

In summary, we have calculated the two-magnon Ra- 
man spectra for various spin models proposed for the 
iron-based superconductors. We have shown that a dis- 
tinct two-magnon Raman response can result between 
models that vary in the level of magnetic frustration. 
Complementary to neutron scattering, magnetic Raman 
scattering detects short-wavelength spin fluctuations and 
can serve as an independent measure of the underlying 
spin interactions. Together with recent experiments, our 
results favor spatially anisotropic models, implying that 
the amount of magnetic frustration in the pnictides and 
chalcogenides is small. 
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FIG. 4. (Color online) S^^{q,oj) [blue peaks] from ED and 
the magnon dispersion spectra [red curves] from spin wave 
theory. Here the momentum points are (in the single iron 
BZ): F = (0, 0), X = (tt, 0), and M = (tt, tt). 
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